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Magnetic  Field,  According  to  Model  3  at  a 
Distance  of  100  km,  1  usec  Rise  Time,  and 
25  usec  Fall  Time  for  the  Ground  and  at 
Five  Altitudes  . 


5.69  Computer  Solutions  to  Equation  (40),  the 
Magnetic  Field,  According  to  Model  3  at  a 
Distance  of  1  km,  2  usec  Rise  Time,  and 
50  usec  Fall  Time  for  the  Ground  and  at 
Five  Altitudes  . 


5.70  Computer  Solutions  to  Equation  (40),  the 
Magnetic  Field,  According  to  Model  3  at  a 
Distance  of  5  km,  2  usec  Rise  Time,  and 
50  usec  Fall  Time  for  the  Ground  and  at 
Five  Altitudes  . 


5.71  Computer  Solutions  to  Equation  (40),  the 
Magnetic  Field,  According  to  Model  3  at  a 
Distance  of  10  km,  2  usec  Rise  Time,  and 
50  usec  Fall  Time  for  the  Ground  and  at 
Five  Altitudes  . 


5.72  Computer  Solutions  to  Equation  (40),  the 
Magnetic  Field,  According  to  Model  3  at  a 
Distance  of  100  km,  2  usec  Rise  Time,  and 
50  usec  Fall  Time  for  the  Ground  and  at 
Five  Altitudes  . 


5.73  Comparison  of  the  Electric  Field  Predicted 
by  the  Three  Models  and  Lin's  Modified 
Model  (Master  et  al. ,  1981)  with  Experimental 
Data  Taken  on  the  Ground  at  a  Distance  of 
10  kilometers.  Estimated  Rise  Time  of 
2  usec  . 
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Figure 

5.74 


5.75 


Comparison  of  the  Electric  Field  Predicted 
by  the  Three  Models  and  Lin's  Modified  Model 
(Master  et  al.,  1981)  with  Experimental  Data 
Taken  at  1500  Feet  at  a  Distance  of  10 
kilometers.  Estimated  Rise  Time  of  2  psec  . 

Comparison  of  the  Electric  Field  Predicted 
by  the  Three  Models  and  Lin's  Modified  Model 
(Master  et  al.,  1981)  with  Experimental  Data 
Taken  at  8300  Feet  at  a  Distance  of  10 
kilometers.  Estimated  Rise  Time  of  2  psec  . 
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Abstract 


A  comprehensive  analysis  and  comparison  of  lightning 
return  stroke  models  is  presented.  A  brief  description  of 
the  lightning  process  and  deviations  of  the  equations  for 
calculating  the  electric  and  magnetic  fields  at  points  on 
the  ground  and  at  points  at  altitude  for  an  arbitrary  cur¬ 
rent  waveform  are  also  presented.  The  models  analysed  are 
the  Bruce-Golde  model,  Uman ' s  transmission  line  model  (the 
breakdown  pulse  current  not  attenuated),  and  Uman ' s  model 
using  an  attenuating  breakdown  pulse  current.  Plots  of  the 
electric  and  magnetic  fields  that  these  three  models  pre¬ 
dict  for  various  distances  and  altitudes  and  different 
parameters  (rise  and  fall  time  of  the  current)  are  included. 
The  above  three  models  along  with  a  recently  modified 
version  of  Lin's  model  (where  the  breakdown  pulse  is  attenu¬ 
ated)  are  then  compared  with  recently  acquired  lightning 
return  stroke  electromagnetic  data.  The  data  was  obtained 
simultaneously  on  the  ground  and  at  various  altitudes  in  a 
WC- 130  aircraft  in  South  Florida. 


AN  ANALYSIS  AND  COMPARISON  OF  LIGHTNING 
RETURN  STROKE  MODELS  AT  ALTITUDE 

I .  Introduction 

A  study  of  how  much  energy  from  a  lightning-sourced 
electromagnet  (EM)  field  couples  into  and  affects  modern 
aircraft  is  necessary  to  insure  continued  flight  safety. 
Modern  aircraft  are  more  susceptible  to  lightning  phenomena 
because  of  tne  increased  use  of  smaller  (and  therefore 
faster,  lower  capacitance,  and  lower  voltage)  electronics 
and  new  carbon  composite  structural  materials  in  their 
design  (Baum,  1980a;  Master,  1981). 

Due  to  their  higher  operating  speeds  and  lower 
capacities,  modern  electronic  devices  are  more  responsive 
or  "resonate"  at  higher  frequencies  and  can  be  destroyed 
at  lower  voltages  than  their  predecessors. 

The  use  of  composite  is  highly  desirable  because 
of  its  lower  cost,  lighter  weight  (and  therefore  lower  fuel 
costs),  and  its  lower  radar  return  signature  (e.g.,  stealth 
aircraft)  as  compared  to  the  more  conventional  use  of 
aluminum.  The  lower  radar  return  these  composite  materials 
offer  is  due  to  their  ability  to  absorb  electromagnetic 
energy  as  opposed  to  reflecting  it  as  aluminum  does.  This, 


however,  also  results  in  a  very  major  disadvantage  in 
using  composite  materials,  that  of  allowing  any  undesirable 
external  EM  energy  into  the  aircraft.  Therefore,  these 
composite  materials  provide  less  protection  or  "shielding" 
of  the  aircraft  interior  to  external  EM  radiation.  This 
means  that  more  energy  from  the  high-frequency,  high- 
voltage  lightning  can  transfer  or  "couple"  into  the 
interior  of  an  aircraft  and  jeopardize  flight  safety  by 
destroying  electronic  components  or  causing  data  and  soft¬ 
ware  errors  in  the  on-board  computers  (Baum,  1980) . 

In  order  to  be  able  to  accurately  calculate  the 
amount  of  energy  due  to  lightning  that  couples  into  an 
aircraft,  one  must  first  have  an  accurate  model  of  the 
lightning  process.  Lightning  has  been  studied  and  modeled 
for  many  years;  however,  it  has  always  been  done  from  a 
surface  of  the  Earth  standpoint.  Not  until  recently 
(Master,  1981;  Pitts  and  Thomas,  1980;  Baum,  1980b)  have 
researchers  been  studying  lightning  at  altitudes  and  try¬ 
ing  to  write  equations  to  predict  the  EM  fields  at  those 
altitudes.  The  only  such  modeling  of  the  EM  fields  pub¬ 
lished  at  this  time  is  one  done  by  Uman,  et  al.  which  uses 
a  slightly  modified  version  of  the  most  recent  model  of 
ground  data  put  forth  by  Lin,  et  al .  (Master,  1981). 
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Problem 


The  problem  of  this  thesis  is  to  investigate  the 
results  of  applying  other  models  for  analysis  and  compari¬ 
son  with  data  taken  at  altitudes  and  determining  how  well 
those  previous  models  predict  the  data.  Also,  one  new 
model  which  utilizes  the  modification  Uman  has  made  to 
Lin's  model  (Master,  1981)  is  investigated. 

Scope 

This  paper  briefly  describes  the  lightning  process, 
presents  the  theory  for  calculating  the  electric  and  mag¬ 
netic  fields,  presents  two  past  models  of  lightning  and 
one  new  model,  presents  the  results  of  those  three  models 
(predicted  values) ,  and  then  compares  the  three  models  with 
actual  data. 

The  three  models  are: 

1.  The  Bruce-Golde  model 

2.  Uman ' s  transmission  line  model 

3.  Uman ' s  transmission  line  model  with  a  height- 
dependent  attenuation 

Each  of  these  models  has  advantages  and  disadvan¬ 
tages  for  use  over  the  other  models. 

Presentation 

Chapter  II  gives  a  brief  description  of  how  the 
lightning  process  occurs  and  defines  pertinent  terminology. 
Chapter  III  derives  equations  for  calculating  the  electric 
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and  magnetic  fields  at  points  on  the  surface  of  the  Earth 
Then  the  more  general  equations  for  calculating  the  EM 
fields  at  altitude  are  presented.  Chapter  IV  discusses 
the  advantages  and  disadvantages  of  each  of  the  three 
models  relative  to  each  other  and  relative  to  Master's 
recent  study  (Master,  1981)  which  uses  a  slightly  modifie 
form  of  Lin's  model.  Chapter  V  presents  the  graphical 
results  of  the  three  models  and  a  discussion  of  how  each 
of  the  models  varies  as  a  function  of  distance,  height, 
and  lightning  parameters  and  compares  each  with  actual 
data.  Chapter  VI  discusses  the  results  of  this  research 
and  makes  some  conclusions  and  recommendations. 
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II .  Background 


Lightning  Process  (Uman,  1969) 

Lightning  is  defined  as  a  "transient,  high-current 
electric  discharge  whose  path  length  is  generally  measured 
in  kilometers  [Uman,  1969:1]."  The  discharge  occurs  when 
one  region  of  space  contains  enough  charge  or  voltage  poten¬ 
tial  relative  to  another  region.  When  the  charge  differen¬ 
tial  is  sufficient  to  cause  electrical  breakdown  of  the 
air  between  the  two  regions,  the  discharge  or  lightning 
occurs . 

The  build-up  of  electrical  charge  occurs  most  often 
in  thunderclouds;  although,  lightning  also  can  be  produced 
by  snowstorms,  sandstorms,  and  the  clouds  over  erupting 
volcanoes.  Sufficient  charge  regions  to  produce  lightning 
have  been  reported  to  exist  even  in  clear  air  (Uman,  1969:1)  . 

As  shown  in  Figure  2.1,  the  major  region  of  posi¬ 
tive  charge  build-up  occurs  at  the  top  of  the  cloud  (P) . 
Moving  down  from  the  top  of  the  cloud,  it  becomes  less 
densely  positive,  eventually  becoming  weakly  negative  and 
then  becoming  strongly  negative  (N) .  At  the  very  bottom 
tip  of  the  cloud,  there  is  a  weak  positive  region  (p) . 

The  lightning  flash  can  be  of  several  types.  The 
most  common  is  called  intracloud  lightning.  As  the  name 
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Fig.  2.1.  Typical  Charge  Distribution  in  a  South 
an  Thundercloud  According  to  Malan  (1952,  1963). 
ric  Field  Measurements  Show  Typically  P=+40  coul , 

N = - 4  0  coul,  and  p=+10  coul  (Uman,  1969  :  3). 


Fig.  2.2.  (a)  The  Luminous  Features  of  a  Ligntning 

Flash  as  They  Would  be  Recorded  by  a  Camera  with  a  Fixed 
Lens  and  Moving  Film.  Time  Scale  Increases  to  the  Right; 
Time  Scale  is  Distorted  for  Clarity.  (b)  The  Same  Lightning 
Flash  as  it  Would  be  Recorded  by  a  Camera  with  Stationary 
Film  (Uman,  1969:6) . 


indicates,  intracloud  lightning  takes  place  within  the 
cloud  itself.  Another  category  of  lightning  is  called 
cloud-to-cloud  lightning.  Other  types  of  lightning  are 
heat,  sheet,  ribbon,  ball,  rocket,  and  bead  lightning 
(Uman,  196  9:  13)  . 

The  type  of  lightning  one  usually  thinks  of  at 
the  mention  of  lightning  and  is  the  most  commonly  studied 
and  written  about  is  called  cloud-to-ground  (CG)  lightning 
sometimes  referred  to  as  streaked  or  forked  lightning. 

This  is  the  type  of  lightning  that  this  thesis  will  inves¬ 
tigate  . 

The  CG  lightning  flash  is  made  up  of  individual 
discharges  called  strokes.  There  is  always  at  least  one 
stroke  per  flash,  typically  there  are  three  or  four.  As 
many  as  26  strokes  have  been  reported  in  one  flash  (Uman, 
1969:10).  Each  stroke,  in  turn,  is  made  up  of  two  parts: 
the  leader  and  the  return  stroke.  The  leader  generally 
propagates  from  the  cloud  to  the  ground  and  is  usually 
negatively  charged.  Relative  to  the  return  stroke,  the 
leader  transfers  less  charge. 

The  first  leader  in  a  flash  is  called  the  stepped 
leader.  The  reason  for  this  terminology  is  that  moving 
film,  fixed  lens  cameras  (e.g.,  a  boy's  camera)  show  that 
this  leader  propagates  in  segments  or  steps  (typically  50 
meter  lengths)  with  pauses  of  approximately  50  usee 
between  steps  (Uman,  1969:5).  A  stepped  leader  would  show 


up  on  a  streak  film  as  a  luminous  step,  a  dark  period,  a 
new  luminous  step  with  a  faint  channel  above  the  new  step 
to  the  cloud,  a  dark  period,  a  new  step  and  faint  channel, 
a  dark  period,  and  so  forth.  The  stepped  leader  is  usually 
strongly  branched  downward  and  is  thought  to  originate  due 
to  electrical  breakdown  between  the  N  and  p  regions  of 
the  cloud  (Uman,  1969:5)  shown  in  Figure  2.1. 

Once  the  stepped  leader  propagates  to  within  a  few 
hundred  meters  of  the  Earth  (ground  or  zero  potential) , 
electrical  breakdown  between  the  high  potential  (usually 
negative  (Uman,  1969:10))  of  the  leader  and  the  ground 
occurs.  When  the  upward  moving  ground  potential  contacts 
with  the  leader,  a  short-circuit-like  situation  exists. 

At  this  point,  an  upward  moving  wavefront  of  ground  poten¬ 
tial  starts  (the  return  stroke).  The  wavefront  travels  up 
the  same  channel  that  the  leader  came  down  (including  any 
"branches"),  eventually  grounding-out  all  negative  charges 
that  were  in  the  leader.  The  return  stroke  is  the  main 
discharge  in  a  flash  and  carries  the  larger  currents. 

The  first  return  stroke  in  a  flash  (i.e.,  the 
return  stroke  following  the  stepped  leader)  will  travel 
up  the  channel  at  a  constant  velocity  until  it  reaches  the 
first  branch  point.  At  this  point,  the  channel  luminosity 
increases  (called  a  branch  point) ,  then  the  return  stroke 
will  continue  at  a  slightly  slower  velocity  until  the  next 
branch  point,  and  so  on.  An  increase  in  the  luminosity 
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of  the  entire  cloud- to-ground  channel  is  known  as  an  M 
component . 

If  the  return  stroke  channel  is  allowed  to  decay 
for  more  than  approximately  100  msec  after  it  reaches  the 
cloud,  the  lightning  flash  is  ended.  If,  however,  more 
negative  charge  is  made  available  to  the  channel  (appar¬ 
ently  from  higher  in  the  N  region  of  the  cloud)  and  the 
channel  is  still  sufficiently  ionized  (typically  less  than 
100  msec  after  the  return  stroke  ends) ,  then  a  secondary 
leader  can  start.  This  new  leader  is  called  a  dart  leader. 
It  will  follow  the  same  main  channel  that  the  stepped 
leader  and  return  stroke  used  (i.e.,  no  "branching"). 

The  dart  leader  appears  to  propagate  an  order  of  magnitude 
faster  than  the  stepped  leader,  but  carries  less  charge. 

Unlike  the  stepped  leader,  the  dart  leader  will 
tend  to  get  closer  to  the  ground  because  the  channel  is 
still  ionized.  Then  a  new  return  stroke  occurs;  however, 
these  secondary  return  strokes  travel  at  a  constant 
velocity,  presumably  due  to  the  lack  of  "branching"  in  the 
dart  leader.  Also,  these  return  strokes  tend  to  have 
higher  velocities  and  transfer  less  charge  than  initial 
return  strokes.  This  leader-return-stroke  sequence  will 
repeat  until  the  charges  in  the  cloud  are  sufficiently 
depleted . 

Strokes  carrying  positive  charge  to  Earth  have  been 
observed  via  direct  measurements  during  discharges  to 
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instrumented  towers  (liman,  1969:10)  in  addition  to  upward 
propagating  leaders  carrying  both  negative  and  positive 
charges  with  no  apparent  return  stroke  (Uman,  1969:11). 
These  types  of  lighting  flashes  will  not  be  investigated 
in  this  thesis. 

Table  2.1  from  the  book  Lightning  (Uman,  1969:4) 
shows  some  of  the  pertinent  parameters  for  the  stepped 
leader,  the  dart  leader,  the  return  stroke,  and  the  light¬ 
ning  flash. 

Assumption 

As  stated  earlier,  this  thesis  is  concerned  with 
investigating  cloud- to-grcund  lightning  flashes;  specifi¬ 
cally,  we  will  calculate  the  EM  fields  at  various  points 
in  space  under  various  conditions  for  each  of  four  models. 
These  calculated  values  will  be  compared  to  experimentally 
obtained  values.  From  this  comparison  a  qualitative  judge¬ 
ment  of  the  "goodness"  of  each  model  performs.  Since  the 
EM  fields  are  of  concern  here,  only  the  return  stroke  is 
of  interest  due  to  its  much  larger  peak  current  relative 
to  the  leader  (thereby  producing  larger  EM  field  values). 

Other  assumptions  used  in  this  study  are  the  fol¬ 
lowing  : 

-9 

1.  Free  space  permitivity  (e^  =  1/36ttx10  )  and 

permeability  (uq  =  4 ::  x  10  ^)  . 


10 


4 


IN 

GJ 

-4 

§ 

£4 


</3 

H 

Z 

M 

Z 

o 

CL, 

a 

o 

u 


cn 

<C 

Cm 

C3 

Z 

M 

Z 

JH 

»-r4 

<3 

H 

4J 


s  ^ 

1  g> 

c  ? 

.5  M 

d  Q) 

2  W 

id  n 
44 

S  O 


03 

U 

•pH 

a 

>, 

H 


C 

c,  03 


O 

O 

IN 


LA 

IN 


O  O 
LD  U") 


O 

CO 


i4 

0) 

03 

<d 

O 

4) 

TJ 

03 

Q. 

Q, 

03 

43 

in 

i 


CL. 

03 

43 

03 

04 

0 

x: 

4-3 

on 

c 

03 

X) 

I 


X 

\D 

IN 


lO 

O 


X 

m 


o 

IN 


in 

o 


X 

o 

H 


04 

U  O 
03 

in  c 
3-  0 


03 

CU 

03 

43 

in 

c 

03 

0) 

? 

43 

03 

-Q 


id 

C3 

01 

Cu 

o 

3-t 


id 

> 

u 

03 

4*3 

c 

•H 

0) 

g 

•H 


a,  w 


44 

o 

>1 

43 

•pH 

o 

o 

rH 

03 

> 

03 

O' 

id 

14 

03 

> 

< 


3 

0 

u 


03 

c 

c 

id 

£. 

u 

c 

o 


u 

03 

x> 

id  -H 


03 


03 


X5 
03 

Qi 

a  o; 


03 

o 

04 

0) 

T> 


0) 

43 

03 


03 

V4 

id 

JO 

u 


r~ 

o 

r— H 

X 


vo 


VO 

O 


X 

o 

fM 


LO 

o 


X 

o 


IM 

O 


O 

03 

133 


3 

0 

O 


03 

C 

c 

id 

30 

u 


43 
id 
tr> 
id  dd 


c 

o 


CL, 

0 


03 

43 


14 

03 

T3 

id 

03 

►4 

43 

U 

10 

Q 


00 

o 


o  o  o  o 

X  oo  m  <-t  la 

A  H  IN 


i— 

o 

rH 


VO 


r- 

o 


X 

o 

IN 


IN 


O 

IN 

I 

O 


o 

•S’ 


O 

03 

0) 


c 

o 


id 

cr> 

id 

a. 

o 

o 


o 

03 

0} 

3- 

id 

x 


03 

0) 

id 

03 

14 

U 

c: 


44 

o 


03 

OS 


u 

03 

03 

OS- 


43 

a 

0) 

u 

u 

3 

U 


•4 

U. 

Qa 

03 

03 

X 

o 

U4 

43 

id 

4H 

Oi 

0 

id 

03 

o 

Q3 

03 

14 

Q, 

>1 

TS 

X 

>1 

aJ 

0 

43 

43 

0 

•H 

03 

14 

*r4 

C 

43 

O 

03 

43 

u 

03 

o 

14 

U3 

0 

14 

0 

rH 

id 

rH 

14 

g 

Q) 

x; 

C 

0) 

3 

■H 

> 

u 

14 

> 

U 

H 

o 

03 

03 

PL 


43 

c 

03 

14 

14 

3 

O 

X 

id 

03 

□4 


id 

X  <43 


43 

C 

a) 

14 

14 

3 

U 

X 

id 

03 

a. 


<d 

jd 

o 

43 

03 

g 

•H 

£4 


11 


•  ,*/»  .  *»  > J 


o 

01 

M 

a 

1 

G 

- 

1 

G 

-P 

P 

O 

C 

O 

(0 

01 

Cn 

01 

p 

CP 

c 

G 

-X 

p 

o 

c 

•H 

rP 

U1 

0 

3 

01 

•r-t 

03 

3 

03 

p 

o 

cn 

Xi 

3 

0 

i — t 

p 

3 

rH 

O 

4-t 

cn 

Cn 

i — 1 

i — t 

0 

G 

G 

U 

3 

X 

V 

P 

c 

■P 

co 

C 

O 

0) 

P 

0J 

01 

3 

G 

-H 

u 

G 

Cb 

01 

C 

r— t 

03 

01 

e 

5 

-P 

4-1 

03 

K 

0) 

1-1 

* 

CO 

-P 

P 

0) 

P 

P 

P 

01 

01 

G 

4-1 

P 

G 

l-i 

3 

V 

a: 

XI 

0 

0 

P 

0) 

01 

U 

x; 

0 

o 

0) 

P 

4-1 

p 

p 

rp 

c 

4-1 

P 

to 

cp 

Cn 

p 

3 

4-t 

0 

CO 

3 

c 

c 

C 

s: 

cn 

> 

0 

■H 

c 

O 

05 

■H 

01 

cn 

P 

p 

03 

P 

3 

rp 

(T3 

44 

0) 

01 

3 

P 

Cn 

P 

G 

f- i 

0 

-p 

u 

P 

p 

C 

•H 

i — 1 

Cn 

G 

c 

3 

•P 

0) 

+J 

01 

p 

*P 

01 

03 

01 

3 

cn  c 

c 

CP 

01 

01 

CP  G 

p 

0 

c 

c 

-9 

aj 

XI 

0) 

P 

■H 

«? 

o 

ns 

-r-l 

a 

G 

s 

05 

4-> 

G 

G 

G 

3 

■M 

-H 

G 

D 

CJ 

P 

2 

e-< 

Eh 

O 

-G 

1 

1 

O' 

| 

t 

1 

1 

•  H 

j 

I 


12 


2.  The  Earth  is  a  flat,  perfectly  conducting, 
infinite  ground  plane. 

3.  Constant  velocity  of  propagation. 

4.  Only  secondary  return  strokes  are  investigated. 

5.  The  channel  is  straight  and  vertical. 

6 .  The  channel  will  be  modeled  as  a  transmission 

1  ine  . 

The  third  and  fourth  assumptions  are  interrelated 
in  that  a  constant  velocity  is  only  a  valid  assumption  on 
analyzing  secondary  return  strokes.  The  choice  of  secon¬ 
dary  return  strokes  prevents  the  difficulty  introduced 
when  trying  to  functionally  describe  the  velocity  of  a 
return  stroke  that  changes  at  every  branch  point  as  is  the 
case  for  first  return  strokes. 

Another  reason  for  studying  only  secondary  return 
strokes  comes  from  the  lack  of  definitive  knowledge  of  the 
EM  fields  due  to  the  stepped  leader/ initial  return  stroke 
process.  Measurements  of  the  first  stroke  show  that  the 
fields  have  a  slower  initial  rise  time  than  secondary 
strokes.  Many  researchers  believe  the  slower  rise  time  of 
first  strokes  is  due  to  the  interaction  of  the  slow,  down¬ 
ward  propagating  stepped  leader  and  the  simultaneous  upward 
moving  ground  leader  before  it  contacts  the  stepped  leader. 
Whereas  dart  leader,  traveling  an  order  of  magnitude  faster 
than  the  stepped  leader,  down  an  already  ionized  channel, 
tend  to  get  closer  to  the  ground  before  the  return  stroke 


begins,  resulting  in  a  faster  initial  rise  time  of  the  EM 
fields.  All  this,  however,  is  not  clearly  understood  and 
therefore  most  return  stroke  studies  deal  only  with  secon¬ 
dary  return  strokes  as  will  be  the  case  in  this  thesis. 

The  channel  will  be  assumed  to  be  straight  and 
vertical.  Some  studies  show  that  path  tortuosity  does 
not  affect  the  induced  EM  fields  (Pearlman,  1979),  whereas 
other  studies  show  that  path  tortuosity  does  have  an 
effect  (Levine,  1979).  The  matter  of  path  tortuosity  is 
as  yet  in  dispute  and  this  thesis  is  not  concerned  with 
resolving  the  issue;  therefore,  because  it  is  easier  to 
model  the  lightning  channel  as  straight  and  vertical,  it 
will  be  so  modeled  in  this  thesis  as  it  is  in  most  other 
studies  not  concerned  with  investigating  path  tortuosity 
effects . 

Lastly,  the  channel  will  also  be  modeled  as  a 
transmission  line  or  antenna  of  finite  length.  Initially, 
in  the  first  two  models,  the  line  will  be  lossless.  For 
the  second  two  models  a  more  realistic  loss  factor  (attenu¬ 
ation)  is  introduced.  Throughout  this  thesis  the  trans¬ 
mission  line  will  have  matched  impedance  loads  at  both 
ends.  Naturally,  the  lightning  channel,  as  can  any  trans¬ 
mission  line,  could  be  modeled  as  a  lumped  parameter  cir¬ 
cuit  (e.g.,  an  LRC  network)  as  is  done  in  a  study  by  Little 
(Little,  1978) .  This  thesis  will  use  more  conventional 
transmission  line  analysis  methods. 
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Ill .  Theory 


EM  on  the  Ground 

Since  it  is  assumed  to  be  straight,  the  lightning 
channel  can  be  looked  upon  as  a  vertical  antenna  with  a 
finite  length.  The  analysis  is  done  using  the  same  methods 
as  used  by  Uman  (Uman,  1975).  Figure  3.1  shows  the 
lightning  model  and  the  parameters  of  concern  for  deter¬ 
mining  E  and  B  at  point  P  on  the  ground. 


Beginning  with  the 

time-dependent  Maxwell's  equa- 

tion 

in  free-space. 

V  •  E  =  p/eQ 

(1) 

V  •  B  =  0 

(2) 

V  x  E  =  -3B/3t 

(3) 

V  X  B  =  U0J  +  If 

c 

(4) 

and 

their  general  solutions 

in  terms  of  retarded  scalar 

and 

vector  potentials. 

E=  -7* -|f 

(5) 

B  =  V  x  A 

(6) 
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where 


4>  (r,  t) 


4^^'SPlr''t‘R/c,dv' 
scalar  potential 


A  (r,  t) 


y  r 

4 ¥./v^  Jfc't  t-R/c)  dv" 
vector  potential 


V 


A  + 


li 

at 


o 


(7) 


(8) 


(9) 


The  scalar  and  vector  potentials  are  based  on  the  geometry 
shown  in  Figure  3.2. 

The  problem  is  to  analyze  a  straight  vertical 
antenna  of  length  H  that  lies  perpendicular  to  a  perfectly 
conducting  ground  plane.  Boundary  conditions  dictate  the 
existence  of  an  image  antenna  of  equal  length  below  the 
ground  plane  as  shown  in  Figure  3.2.  The  antenna  cross- 
section  is  much  less  than  the  radiation  wavelength.  The 
current,  i(z,t),  is  a  continuous  function  and  is  zero 
everywhere  for  t£0 . 

.  —  A 

Since  i(z,t)  lies  along  the  z-axis,  r'  =  zaz 

and 


J ( r t-R/c) dv"  =  i ( z , t-R/c) dzaz 


(10 


To  derive  the  differential  magnetic  field,  dBR, 
to  an  infinitesimal  current  element  on  the  real  antenna. 


du' 


Fig.  3.2.  The  Geometry  for  General  Solutions  of 
the  Time-Dependent  Maxwell  Equations  (Uman,  et  al . ,  1969 


start  with  Maxwell's  equations.  Equation  (6)  becomes 


dBR  =  7  X  dAR 

To  find  dA  ,  substitute  (10)  into  (8)  to  get 

i\ 

—  —  /*H  1 

Ar  ( r  / =4 ~  Jq  r  Mz,t-R/c)azdz 


which  means 


dAR  =  47  |  i(z,t-R/c)azdz 

To  convert  Equation  (13)  from  cylindrical  to  spherical 
coordinates,  let 

a  =  a  cos  6  -  a,  sin  9. 
z  r  6 

Equation  (13)  becomes 

dA_  =  tt  4  i(Zft-R/c)  [a  cos  6-  a,  sin  9]  dz 

K  4  11  IT  t? 

=  dA  a  +  dA~a„  +  dA  a, 

r  r  u  O  9  4> 

Substitute  Equation  (14)  into  Equation  (11)  to  get 


r  R  sin 


(dAb  sint,)  ~  <dA^>  1 


1  1  3  3 

afi  4  L  1  ■  M  yr  (dA_)  -  --y  ( RdA*  )  ] 

u  R  sm  b  d $■  R'  or  v 


o  3 

+  a,  [-^  ( RdA^ )  -  ^  ( dAp )  ] 

y  o  R  v,  do  R 


By  comparing  Equations  (14)  and  (15),  it  is  seen  that 
cLV  =  0.  Also,  dAra  has  no  dependence  on  d;  therefore 
derivatives  with  respect  to  1  are  zero  (i.e.,  S/Dd  = 


Equation  (15)  becomes 


_  r1 

dIR  =  I  %  ["R  4?  I  i  ( z  ,  t-R/c)  sin  £  dz] 


~  jk  [4^  I  i(z, t-R/c)  cos  6  dz]} 


Or  1  .  -  9  .  ,  „  „  .  . 

=  a,  -7—  l*rSin  t  t—  i(z, t-R/c) 

i  4  TT  R  oR 


+  —  sin  5  i  (  z  ,  t-R/c )}  dz 


Using  the  identity 


—  i(z, t-R/c)  = 

-•R 


Equation  (17)  becomes 


JR  4tt 


i(z, t-R/c) 


i ( z , t-R/c) 


i ( z , t-R/c) ]  dz 


■»-v  *  -V  m  C"  V'T-'T  f 


•t  ■!  ■  i  ■  i  ■■i  1 


■  V.  y  r* 


c 


re 


As  mentioned  earlier,  boundary  conditions  require 
an  imaginary  differential  current  element  with  the  same 
magnitude,  direction,  and  distance  below  the  ground  plane 
as  the  real  current  element  above  the  plane.  Now  proceed 
to  calculate  the  differential  magnetic  field,  dB,.,  due  to 
infinitesimal  current  element  as  was  cone  or  L-he  real 
current  element  above  to  find 


d3T  =  a.-rm  sin  &  [-K  i(z,t-R/c) 
1  -4"  R2 


+  ~  i(z,t-R/c)]  d: 

CK  c  l. 


=  dB 


R 

The  total  differential  magnetic  field  is  the  sum 
of  the  fields  due  to  the  real  and  imaginary  currents. 
Thus 


d3  =  dB0  +  dB-j. 


sin  3  i(z,t-R/c) 

y  TF  y-\  ~ 

rl 


-  cR  7E  i(Z't"R/C)]  dZ 
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Find  the  electric  field  due  to  the  infinitesimal 
real  current  by  starting  with  Equation  (5) 


21 


is  needed. 


Get  it  by  solving  Equation  (9) 


First,  dtp 


to: 


•  A 


+ 


1_  £0 
2  C't 


0 


c 


2  n 
-c  V 


A  •’  t 


(21 


Since  Equation  (9)  involves  the  derivative  of  q>  with 
respect  to  time,  integrate  Equation  (21)  over  time  to  undo 
the  derivative. 


*  ( R / 1) 


A  dt 


The  differential  scalar  potential  due  to  the  real  current 
element  being 


d; 


R 


( R  t 


t) 


(22 


Taking  the  divergence  of  Equation  (14)  results  in 

7  •  dAR  =  ~  cosi[-^  [Ri  ( z ,  t-R/c)  ] 

R“ 


-  2i(z, t-R/c) ; 

Substituting  this  into  Equation  (22)  and  noting  that 
r  c“  =  1/e  yields 


i 


a:,(R,') 


cos  0  [—■  / 

'0  R“  -/0 


i ( z , t -R/c)  do 


i ( z , t-R/c) 


.ho  negative  gradient  of  Equation  (23)  to  get 


1  "> 

d;„  =  a  — -  cos  a  [— — =•  i(z, t-R/c) 


R"  J0 


i  (  z  ,  t -R/c)  d~ 


(  z , t-R/c)  ]  dz 


a.  -.-A-  4  sin  5  l~  /  i  (  z ,  t-R/c)  dt 
'  4  un  R  -'0 


+  i  (  z  ,  t-R/c)  ]  dz 


Using  Equation  (14),  take  the  derivative  with  respect  to 
tine  and  cet 


-k  dAP.  "  ar  rr  !lcos;  -k  i  V, t-R/c)  ]dz 


a-  A-  [4  sinC  A-  i  ( z  ,  t-R/c) ] dz 
4  .  k  o  t 


.nally,  to  find  the  differential  electric  field  due  to 
io  real  current  element,  add  Equations  (24)  and  (25) 
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i ( z , T -R/c) dt 


+  — i  ( z , t-R/c) ] a 
cR" 

1 

+  sin  6  [ — y  j  i ( z , t-R/c) dx 
R  JQ 


+ 


CR 


i ( z , t-R/c)  + 


c2R 


d  t 


i ( z , t-R/c) ] a , } az  (26) 


The  electric  field  due  to  an  infinitesimal  imaginary- 
current  element  can  be  found  in  a  similar  manner  with,  as 
shown  in  Figure  3.1,  v  replaced  by  it  -  G .  The  result  is 

dE  =  .  {  -  cos  9  [~  f  i(z,x-R/c)dx 

0  R  J0 


+  — ^  i(z, t-R/c) ] a 
cR  rl 


1  f 

+  sin  6  [— y  I  i ( z, t-R/c) d^ 
P.  *'0 


+  — i(z,  t-R/c)  +  - - r-  i(z, t-R/c) ]a,  }dz 

cR~  c  R  I 

(27) 

Before  adding  Equations  (26)  and  (27)  to  get  the 
total  electric  field,  the  unit  vectors  a^ ,  a.,  a„  ,  and 

i\  C  IN.  ^ 

ap  need  to  be  converted  to  common  terms.  The  four  unit 

vectors  can  be  broken  up  into  two  components:  a  vertical 

component  (a  )  and  a  horizontal  component  (a) .  From 
z  u 

Figure  3.1  it  is  seen  that 
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< 


< 


a  - 
r 


■a  a  +  aITaM 
z  z  H  H 


=  -  sin  ( :  -  a  +  cos  (6  -  — )  a._ 
2  z  2  H 


=  cos  t-  a  +  sin  t  aT, 
z  H 


a0  =  ~azaz  "  aHaH 


=  ~  c os  (  ?  -  a  -  s  in  ( 5  -  ^ 5  ar_ 
2  z  ^  h 


=  -  sin?  a  +  cos  r  a,, 
z  H 


a  =  a  a  +  aTIatI 
z  z  H  H 


=  cos  ("-r)  a  +  sin  (  r-9)  a„ 

Z  ri 


-  cos-  a  +  sm  -  a 


H 


a  =  -a  a  +  aTIaTI 
j  z  z  H  H 


=  -  sin  (  ”- i)  a z  +cos('-'i)a 


=  -  sin i  a  -  cos?  aTT 
z  H 


The  total  differential  electric  field  is  the  sum 


ields  due  to  the  real  and  imaginary  currents 


•  '  'V  I  ■*. 1  ■  V  1  I'M"  ■  1  ■  ■ 


k 


i 


Thus,  with  appropriate  substitutions  of  Equations  (28), 
(29) ,  (30) ,  and  (31) 


dE  -  dE  +  dE 
K.  X 


4TTS 


{cose  [ 


-X  f  i(z,T- 
R  J0 


R/c) dT 


+  - j  i ( z , t-R/c) ] (cosSa  +sin5a  ) 

cR  Z 


+  sin  0  [-Xj-  f  i  ( z ,  t-R/c)  d  l  +  -  i(z, t-R/c) 
R  J0  cR~ 


+  - X—  -X-  i ( z , t-R/c) ]  (-sin  6  a^+cos  6  a  ) )  dz 


2d  Ot 
c  R 


~r —  {-COS  fi  [-^r  f  i(z,  T-R/c)  di 
R  -'0 


+  - ^  i  (  z  ,  t-R/c)  ]  ( -  cost  a  +  sint  a.  ) 

,  Z  Z  tl 


cR 


+  sin 


rt  i 

I  i (z, t-R/c) dt  +  — ^  i(z, t-R/c) 
JQ  cR“ 


• 

+  ?  Jt 

c  R 

-  1 

2  :  ~  0 

• 

? 

+  0  i  ( 

cR2 

R3  -'O 


i  (  z ,  t-R/c)  ]  (  -  sint  a  -  cost  au)  )  dz 

Z  n 


2  r  ,  2 

COS 


2  f^~ 

— ^  /  i(z, T-R/c)  dT 

R  ^0 


}  ■ 
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U  s ing  t 


+ 


- ~  i(z,t-R/c)  +  —5—  —  i ( z , t-R/c) ] a  j dz 

cR2  cR  2 


Tr-i—  [^r  (2  cos26  -  sin2  0)  /  i  (  z  ,  t-R/c)  di 

27,C0  R° 


(  2  cos^S  -  sin“  3 )  i  {  z  ,  t-R/c) 

cR2 


— h—  ■—  i(z,t-R/c)]a  dz 
c  R  ct  2 


rig  identities 

o  1  1  1 

2  cos  v  -  sin“£  =  2  x  —(1  +  co  s26)  -  -(l-cos22) 

=  ^  cos  29 

=  ^ ( 1  +  3  cos29  ) 

=  ~(1  +  3 (1-2  sin2  r  )  ) 

=  ^(1  +  3-6  sin2  j) 

=  2  -  3  sin  v 


Substituting  Equation  (33)  into  Equation  (32) 


dE  = 


yrr  (2-3  sin  '-’)  /  i(z,x-R/c)dT 

0  R  -0 


+  — (2-3  sin2  9)i(z,t-R/c) 
cR" 


sin^'3^— •  i(z,t-R/c)]a  dz 
c  “R  J  2 


To  get  the  magnetic  and  electric  fields  due  to 
the  current  along  the  entire  length  of  the  channel,  inte¬ 
grate  Equations  (19)  and  (34),  respectively,  from  zero  to 
maximum  height,  H. 


B,(D,t)  =  f  dB,(E,9,t)dz 

r  J  0 


/***  '•‘0  1 
=  /  grr  Sin  e  [—  i  (z,t- 

«/  n  **  "  o  “ 


R/c) 


+  ^  7E  i<z,t-R/c)]dz 


n  /*  R  i 


0  R 


0  /* j_ 

=  Trn'  /  — ir  sinf  i(z,t-R/c)dz 

2"  Jo  R2 

+  2f/  ^  ^  iU.t-R/cldz 


iJo  rH 


E  (D 
z 


,  t)  =  r  dE  (R,3,t)dz 
-/0  2 


1  1  o 

— —  f^er  (2-3  sin" 


•>  /_t  i  <  ~ ~  ~ 


R/c )  d 


The  first  term  in  Equation  (35)  is  called  the 
induction  or  intermediate  term.  The  second  term  is  called 
the  radiation  or  far-field  term.  The  terms  in  Equation  (36) 
are  called  the  electrostatic,  the  induction,  and  the  radia¬ 
tion  terms. 

Because  of  the  integral  of  the  current,  the  electro¬ 
static  term  should  have  the  longest  lasting  effect  on  the 
electric  field.  It  should  result  in  a  theoretically  con¬ 
stant  value  (actually  it  would  slowly  decay)  in  the  elec- 

3 

trie  field.  Due  to  the  1/R  dependence,  the  electrostatic 

term  has  less  effect  the  further  away  the  calculated  point 

is  from  the  return  stroke. 

The  induction  terms  also  have  less  effect  on  the 

EM  fields,  but  do  not  drop  off  as  fast  as  the  electrostatic 

2 

term  since  the  induction  term  is  dependent  on  1/R  ; 
therefore,  the  induction  term  has  a  longer  effect  on  the 

n  n 


field  (in  terms  of  distance  away  from  the  return  stroke) 
than  does  the  electrostatic  term.  Because  it  is  a  scaled 
value  of  the  current,  the  induction  terms  should  produce 
a  hump  in  the  field. 

The  radiation  terms  have  the  most  effect  on  the 
fields  the  further  from  the  return  stroke  one  gets.  This 
is  due  to  it  having  only  a  1/R  dependence.  Also,  due  to 
its  partial  derivative  with  respect  to  time  dependence,  it 
v.’ill  be  predominate  in  the  first  few  hundred  nanoseconds 
(until  the  peak  of  the  current)  after  the  return  stroke 
starts . 

CM  Everywhere 

Equations  (25)  and  (36)  are  good  only  for  calcu¬ 
lating  the  electric  and  magnetic  fields  at  points  on  the 
ground  plane.  Using  methods  similar  to  deriving  Equations 
(35)  and  (36),  equations  for  determining  the  E  and  B  at 
any  point  in  space  can  be  derived. 

Looking  at  Figure  3.3,  it  is  seen  that  the  dis¬ 
tance  from  the  real  current  element  to  the  field  point, 

R  ,  and  the  distance  from  the  imaginary  current  element 
to  the  field  point,  R  ,  are  no  longer  always  equal.  Also, 
the  inside  angle  between  R  and  the  vertical  is  no  longer 
'•  -  .  Taking  these  two  differences  into  account,  Equation 

(26)  with  Equations  (28)  and  (29)  substituted  and  noting 
that  sin  i  =  D/RR  and  cos  •;  =  (HA-r)/P.R  becomes 


30 


--5 - -  j  i(2,T-R  /C)d- 

rr  ^ 


,  2  (HA-z)  2-D2  .  ,  ,  „  ,  , 

+  - 4 -  i(z,t-R  /c) 

cRr 


D 

— - m  -^r  i(z,t-R  /c)  ]az 

c“RP 


+  3D(Hg--Z-).  f  i(Z,T-R^/C)dT 

Rr3  7° 


-(HVLl  i(2»t-R  /c) 
cRR 


+  ^  i(z,t-RR/c)]ar}dz 

c  rr 


Comparing  Equation  (37)  with  the  equations  found 
in  Master,  et  al.  (Master,  1981),  it  is  seen  that  Equation 
(37)  is  equivalent  with  different  notation. 

The  differential  magnetic  field  for  all  points  in 
space  is  simply  found  by  letting  sinf  =  D/R  in  Equation 
(18)  to  a  e  t 


d3R  =  2j~r  [-^j  i(z,t-RR/c) 
rr 


-At  i ( z , t-R„/c) ] dza . 


which  is  equivalent  to  the  equations  found  in  Master,  et  al . 
(Master,  1981)  w7ith  appropriate  notational  changes. 

The  electric  and  magnetic  fields  due  to  the  differ¬ 
ential  length  image  current  element  is  found  by  substituting 
Rt  for  and  -z  for  z  into  Equations  (37)  and  (38) .  The 
resulting  dE^  and  dE^  are  added  to  get  the  total  differen¬ 
tial  electric  field  dE.  Likewise,  dB  is  obtained  by  adding 
d3D  and  d3T. 

The  E  and  B  are  found  by  integrating  dE  and  dB  from 
zero  to  H,  respectively.  Finally,  it  is  found  that 


( D ,  HA ,  t )  =  [fH  3-D-lH-|  *1  f  i(z,T-R  /c)didz 

^  n  J  c\  K 


■  0  ^  0  R 


R 


L 


-H 


3DUIA+Z  )  rt 

—  /  i(-z  ,t-R  /c)dtdz. 

•'0  1  1 


0  V 
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where  Rri  =  (D“+  (HA-z)  “)  t 
R 

and  R,  =  (D2+ (HA+z) 2)  5 
1 

Just  as  was  seen  in  the  E  and  B  equations  for 
ground  point  calculations,  the  eleceric  field  has  an  electro 
static  term,  an  induction  term,  and  a  radiation  term  in  both 
the  vertical  and  horizontal  directions.  Obviously,  since 
it  is  the  same  equation,  the  magnetic  field  still  has  an 
induction  term  and  a  radiation  term  in  the  c-direction. 

The  characteristics  of  these  terms  as  discussed  in  the  pre¬ 


vious  section  still  hold. 


IV.  Models 

As  was  mentioned  i;  Chapter  I,  this  thesis  investi¬ 
gates  three  models  of  the  current  in  the  lightning  channel. 
They  are: 
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It  consists  of  a  uniform  current  that  propagates  up  the 
channel.  This  uniform  current  has  a  time-dependent  decay 
factor  (in  the  form  of  a  double-exponential)  applied  to  it. 
The  BG  model  can  be  thought  of  as  returning  charge  stored 
in  the  corona  envelope  surrounding  the  lightning  channel  to 
ground  (Golde,  1977:331;  Lin  et  al . ,  1980).  Figure  4.1 
pictorially  displays  this  model. 

An  advantage  of  the  BG  model  is  that  its  simplicity 
provides  easy  calculations  of  the  EM  fields.  A  disadvan¬ 
tage  of  the  BG  model  is  that  it  is  physically  untenable 
because  it  has  a  discontinuity  at  the  leading  edge  of  the 
current.  Also,  the  time-dependent  decay  results  in  the 
current  channel  between  the  leading  edge  and  the  ground 
changing  uniformally  and  equivalently  throughout.  This 
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means  the  channel  must  be  able  to  instantaneously  transfer 
charge  from  any  point  along  the  channel  to  ground. 

Another  disadvantage  to  the  BG  model  is  that  the  velocities 
and  peak  currents  it  predicts  seem  to  be  distant  dependent 
(Lin,  1978:128;  Lin  et  al . ,  1980). 

The  second  model  was  first  proposed  by  Uman  and 
McLain  in  1969.  Uman  and  McLain's  model  is  commonly  known 
as  the  transmission  line  model  because  the  negative  z/v 
factors  in  the  exponentials  cause  the  pulse  to  propagate 
up  the  lightning  channel  as  it  would  up  a  transmission 
line.  (Uman  and  McLain  used  a  non-exponential  pulse  in 
their  model,  but  in  this  paper  the  pulse  will  be  modeled 
as  a  double-exponential.)  The  current  in  this  model  does 
not  represent  a  removal  of  charge  from  the  corona;  rather, 
it  represents  an  impulse-like  current  due  to  the  atmos¬ 
pheric  breakdown  between  the  end  of  the  leader  and  the 
start  of  the  return  stroke.  Figure  4.2  shows  the  current 
representation  in  this  model . 

The  advantage  to  the  transmission  line  model  is 
that  it  is  more  physically  reasonable  than  the  BG  model; 
there  is  not  any  discontinuity  of  the  current  at  the  lead¬ 
ing  edge  nor  are  there  any  instantaneous  charge  transfers 
through  the  channel.  Another  advantage  of  the  transmission 
line  model  is  that  it  provides  easy  analytical  computation 
of  the  EM  fields  relative  to  more  recent  models  (i.e., 

Lin's  model,  1978). 
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Lin  (1978)  showed  that  a  primary  disadvantage  to 
Uman  ana  McLain's  model  is  that  it  requires  the  propagating 
pulse  to  have  a  much  shorter  duration  than  direct  measure¬ 
ments  showed  to  be  the  case  (Lin,  1978:128) .  Another  dis¬ 
advantage  of  the  transmission  line  model  is  that  it  pre¬ 
dicts  peak  currents  generally  greater  than  30  kA  and  often 
as  high  as  50  to  100  kA;  however,  direct  measurement  shows 
the  peak  current  to  be  usually  on  the  order  of  10  kA 
(Lin  et  al.,  1980)  .  Yet  another  disadvantage  of  this  model 
is  that  for  far-fields,  it  predicts  the  occurrence  of  a 
"mirror  image"  effect  that  rarely  occurs  in  the  data 
(Master  et  al . ,  1981). 

The  third  model  is  the  same  as  the  second  model, 
but  the  breakdown  pulse  undergoes  attenuation  as  it  propa¬ 
gates  upward  through  the  channel.  The  reason  for  adding 
this  attenuation  is  twofold:  (1)  In  1980,  Jordan  and  Uman 
showed  that  the  channel  luminosity  (and  hence,  by  implica¬ 
tion,  the  current)  decreased  with  height.  (2)  Allowing 
the  breakdown  pulse  to  attenuate  to  a  negligible  value 
when  it  reaches  the  top  of  the  channel  should  greatly  reduce 
the  undesirable  "mirror  image"  effect  (Master  et  al.,  1981). 

Because  the  pulse  current  is  modeled  as  a  double- 
exponential,  the  current  function  in  the  third  model  becomes 
very  close  to  the  corona  current  that  Lin  uses  in  his  model 
(explained  below) .  The  differences  are  that  in  the  third 
model  in  their  thesis  the  time  constants  in  the  double- 
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exponential  (a  and  £)  are  such  that  a  much  faster  rise  and 
fall  time  occurs  and  that  the  corona  current  in  Lin's  model 
(and  the  BG  model)  does  not  propagate  (no  z/v  factors  in 
the  exponentials) . 

The  advantage  to  this  model  is  that  it  is  some¬ 
what  more  physically  realistic  as  outlined  above;  essen¬ 
tially,  the  transmission  line  is  now  no  longer  lossless. 

As  mentioned  in  Chapter  I,  one  model  has  been 
analytically  examined  at  altitude.  That  model  is  a  modi¬ 
fied  form  of  Lin's  model.  Lin's  model  consists  of  three 
currents : 

1.  breakdown  pulse  current, 

2.  uniform-  current,  and 

3.  corona  current  that  "turns  on"  after  the  peak 
of  the  breakdown  pulse  current  has  passed  (Lin,  1978). 
Originally,  the  breakdown  pulse  current  did  not  undergo  any 
attenuation.  In  this  fourth  model,  however,  the  breakdown 
pulse  current  is  attenuated  as  it  is  in  the  third  model 
(Master  et  al . ,  1981).  Fig.  4.3  pictorially  represents 
Lin ' s  model . 

The  advantage  to  Lin’s  model  is  that  it  has  been 
shown  (Lin,  1978)  to  more  closely  predict  the  EM  fields 
than  the  BG  or  transmission  line  models.  The  disadvantage 
is  that  it  is  much  more  analytically  complicated  to  predict 
the  EM  fields  with  this  model. 
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Slope r 


Looking  at  the  above  models,  one  notices  that 
there  are  generally  three  physical  characteristics  of  the 
lightning  being  taken  into  account  in  the  form  of  currents. 
These  three  currents  are: 

1.  current  due  to  the  atmospheric  breakdown, 
having  an  impulse-like  structure,  I  ; 

2.  constant  uniform  current  due  to  the  leader,  I  ; 


and 


3.  current  due  to  the  charges  stored  in  the  corona 
envelope  going  to  ground,  Ic. 

Also,  the  currents  have  three  types  of  attenuation 
or  decay: 

1.  no  attenuation  or  decay  te.g.,  the  second  model) 

2.  height-dependent  attenuation  [exp(-z/X)] 

(e.g.,  the  pulse  current  in  the  third  and  fourth  models); 
and 

3.  time -dependent  decay  [exp(-at) -exp(-£t)  ] 

(e.g.,  the  corona  current  in  the  first  and  fourth  models). 

The  first  and  second  models  have  been  analysed  and 
compared  only  at  ground  level.  The  third  model  has  never 
before  been  analyzed  or  compared  with  data.  The  fourth 
model  has  been  analyzed  at  altitude  (Master  et  al . ,  1981) 
but  has  yet  to  be  compared  to  any  real  data. 
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V .  Analysis 

Theoretical 

In  the  theoretical  analysis,  computer  solutions  to 
Equations  (39)  and  (40)  for  each  of  the  three  models  were 
found.  Figures  5.1  through  5.36  are  the  electric  field 

predictions  and  the  units  on  the  vertical  scale  for  all  of 

them  is  volts  per  meter.  Figures  5.37  through  5.72  are 
the  magnetic  field  predictions  with  units  of  webers  per 
square  meter  on  the  vertical  scale. 

For  each  model,  the  rise  and  fall  time  for  the 
pulse  or  double  exponential  was  varied.  Each  model  was 
analyzed  for  the  following  rise  and  fall  times: 

1.  500  nsec  rise  time,  10  usee  fall  time 

2.  1  usee  rise  time,  25  usee  fall  time 

3.  2  usee  rise  time,  50  usee  fall  time 

For  each  rise  and  fall  time,  the  distance  was 
varied.  The  distances  selected  were  one  kilometer,  five 
kilometers,  ten  kilometers,  and  one  hundred  kilometers. 

Finally,  at  each  distance,  the  fields  were  calcu¬ 
lated  at  ground  level,  fifteen  hundred  feet,  five  thousand 
feet,  eight  thousand  feet,  fifteen  thousand  feet,  and  thirty 
thousand  feet.  All  six  altitudes  for  each  distance  are  in 
one  figure.  The  top  left  graph  in  each  figure  represents 


ground  level;  the  center  left  graph  represents  fifteen 
hundred  feet;  the  bottom  left  graph  represents  five  thousand 
feet;  the  top  right  graph  represents  eight  thousand  feet; 
the  center  right  graph  represents  fifteen  thousand  feet; 
the  bottom  right  graph  represents  thirty  thousand  feet. 

The  distances  and  altitudes  used  were  chosen  to 
correspond  with  both  experimental  data  (taken  in  August 
1981)  and  with  the  predictions  by  Master  (Master  et  al., 

1981)  . 

Table  5.1  shows  the  breakdown  of  the  electric  and 
magnetic  field  figures  for  the  various  parameters. 

In  any  of  the  plots,  the  radiation  term  predominates 
at  the  start.  This  is  because  the  fast  rise  time  of  the  cur¬ 
rent  yields  a  large  positive  value  for  its  derivative  with 
respect  to  time. 

The  initial  radiation  effect  in  the  first  model  is 
generally  shorter  in  duration  and  magnitude  than  in  the 
second  and  third  models.  This  is  due  to  the  pulse  not 
propagating  in  the  first  model.  Initially,  the  radiation 
term  predominates  because  of  the  fast  rise  time  of  the  cur¬ 
rent;  however,  the  derivative  is  at  its  maximum  at  the 
beginning  and  immediately  starts  dropping  off,  eventually 
going  to  zero.  As  the  radiation  term  drops  off,  the  induc¬ 
tion  term  (and  the  electrostatic  term  in  the  electric 
field)  is  getting  larger.  In  the  first  model,  the  current 
is  represented  by  a  propagating  uniform  current  that 
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changes  amplitude  according  to  a  double  exponential; 
nowhere  along  the  channel  at  any  instant  of  time  can  this 
double  exponential  shape  be  seen  (and  hence,  the  rise 
time) .  In  the  second  and  third  models,  however,  the  double 
exponential  shape  is  always  seen  along  the  channel  (until 
it  reaches  the  top)  since  these  two  models  use  a  propa¬ 
gating  pulse.  Consequently,  integrating  from  the  ground  to 
the  height  of  the  return  stroke  at  any  instant  of  time  would 
not  include  the  fast  rise  time  in  the  first  model  whereas 
it  would  be  included  in  the  second  and  third  models. 

The  difference  between  the  second  and  third  models 
is  that  the  'nitial  peak  or  effect  of  the  radiation  term 
is  slightly  less  in  the  third  model  because  of,  the  height 
dependent  attenuation  (e  '  ') .  In  the  first  few  hundred 
nanoseconds  during  which  the  radiation  term  predominates, 
the  pulse  has  not  propagated  very  far  up  the  channel  and 
consequently  the  attentuation  has  little  although  notice¬ 
able  effect  at  this  point. 

In  the  work  done  by  Master  et  al.  using  the  fourth 
model  (Master  et  al . ,  1981),  it  was  seen  that  attenuating 
the  breakdown  pulse  did  not  much  affect  the  fields  on  the 
ground  at  all  distances  and  the  fields  at  all  attitudes 
beyond  ten  kilometers.  This  was  due  to  the  radiation  term 
predominating  the  initial  portion  of  the  fields  while  the 
pulse  is  near  the  ground  and  not  having  much  attenuation. 
Latter  portions  of  the  fields  are  due  to  the  uniform  and 
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corona  currents  which  did  not  change  from  Lin's  original 
model  (Master  et  al. ,  1981)  . 

In  comparing  the  second  and  third  models,  it  is 
seen  that  there  is  quite  a  difference  between  the  two 
models  at  all  distances  and  all  altitudes.  This  would  be 
due  to  there  being  no  uniform  and  corona  currents  to  pre¬ 
dominate  later  portions  of  the  fields. 

Another  difference  between  the  two  models  is  that 
the  third  model  has  essentially  gotten  rid  of  the  undesired 
"mirror  image"  effect.  This  is  due  to  the  pulse  having 
been  greatly  attenuated  by  the  time  it  reached  the  top  of 
the  channel . 

Looking  at  Figures  5.1  through  5.72,  it  is  seen 
that  the  positive  portions  of  the  electric  field  decrease 
while  the  negative  portions  increase.  This  is  especially 
true  in  the  second  and  third  models.  Eventually,  above  a 
certain  altitude  (higher  with  increasing  distance),  the 
field  becomes  essentially  totally  negative.  This  is  appar¬ 
ently  due  to  more  of  the  channel  lying  below  the  local 
horizon  of  the  observation  point. 

In  the  plots  of  the  electric  and  magnetic  fields 
predicted  by  the  second  and  third  models,  it  is  seen  that 
at  a  given  distance  as  altitude  increases,  the  maximum 
magnitude  increases  while  at  the  same  time  the  field  becomes 
more  and  more  impulsive.  This  could  be  due  to  a  decrease 
in  the  effect  of  the  radiation  term  due  to  the  time 
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Fig.  5.5.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  According  to  Model  1  at  a  Distance  of 
1  km,  1  usee  Rise  Time,  and  25  usee  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Tic.  5.8.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  According  to  Model  1  at  a  Distance  of 
100  km,  1  ;:sec  Rise  Time,  and  25  usee  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fig.  5.11.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  According  to  Model  1  at  a  Distance  of 
10  km,  2  ..sec  Rise  Time,  and  50  usee  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fiq.  5.12.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  According  to  Model  1  at  a  Distance  o 
100  km,  2  w sec  Rise  Time,  and  50  ..sec  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fig.  5.13.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  According  to  Model  2  at  a  Distance  of 
1  km,  500  nsec  Rise  Time,  and  10  ::sec  Tall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fiq.  5.14.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  Accordinq  to  Model  2  at  a  Distance  of 
5  km,  500  nsec  Rise  Time,  and  10  -..sec  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fia.  5.15.  Computer  Solutions  to  Equation  (39), 
he  Electric  Field,  Accordinct  to  Model  2  at  a  Distance  of 
0  km,  500  nsec  Rise  Time,  and  10  ..sec  Fail  Time,  for  the 
round  and  at  Five  Altitudes. 
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Fin.  5.16.  Commuter  Solutions  to  Cauation  (39), 
ctric  T’ield ,  rccordina  to  Model  2  at  a  Distance  of 
500  nsec  Hi  so  Time,  and  10  -..sec  Fall  Time,  for  the 
:i n d  at  r i  ve  2 ■.  1 1  itu d o s . 
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Fig.  5.17.  Computer  Solutions  to  Eauation  (39), 
the  Electric  rield,  Accordina  to  Model  2  at  a  Distance  of 
1  km,  1  .isec  Rise  Time,  and  25  nsec  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fig.  5.18.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  Accordina  to  Model  2  at  a  Distance  of 
5  km,  1  ::rec  Rise  Time,  and  25  usee  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fig.  5.19.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  Accordino  to  Model  2  at  a  Distance  of 
10  km,  1  i.sec  Rise  Time,  and  25  ;;sec  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fiq.  5.20.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  According  to  Model  2  at  a  Distance  o 
100  km,  1  usee  Rise  Time,  and  25  usee  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fig.  5.21.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  According  to  Model  2  at  a  Distance  of 
1  km,  2  rsec  Rise  Time,  and  50  ;:sec  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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5.22.  Computer  Solutions  to  Equation  (39), 
Field,  According  to  Model  2  at  a  Distance  of 
Rise  Time,  and  50  usee  Fall  Time,  for  the 
Five  Altitudes. 
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Fiq.  5.23.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  According  to  Model  2  at  a  Distance  of 
10  km,  2  nsec  Rise  Time,  and  50  nsec  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fig.  5.24.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  According  to  .Model  2  at  a  Distance  of 
100  km,  2  ysec  Rise  Time,  and  50  ,isec  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fig.  5.25.  Computer  Solutions  to  Equation  (39), 
The  Electric  Field,  Accordinq  to  Model  3  at  a  Distance  o 
1  km,  500  nsec  Rise  Time,  and  10  usee  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fig.  5.26.  Computer  Solutions  to  Equation  (39), 
The  Electric  Field,  Accordinq  to  Model  3  at  a  Distance  o 
5  km,  500  nsec  Rise  Time,  and  10  usee  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fiq .  5.27.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  According  to  Model  3  at  a  Distance  of 
10  km,  500  nsec  Rise  Time,  and  10  usee  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fig.  5.28.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  According  to  Model  3  at  a  Distance  of 
100  km,  500  nsec  Rise  Time,  and  10  usee  Fall  Time,  for  th 
Ground  and  at  Five  Altitudes. 
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Fig.  5.29.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  According  to  Model  3  at  a  Distance  of 
1  km,  1  iisec  Rise  Time,  and  25  usee  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fiq.  5.30.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  According  to  Model  3  at  a  Distance  of 
5  km,  1  :;sec  Rise  Time,  and  25  nsec  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fiq.  5.31.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  Accordina  to  Model  3  3t  a  Distance  o 
10  km,  1  ,:sec  Rise  Time,  and  25  ..sec  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fiq.  5.32.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  According  to  Model  3  at  a  Distance  o 
100  km,  1  ..sec  Rise1  Time,  and  25  ;.sec  Fall  Time,  for  the 
Ground  and  at  Five  Altitudes. 
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Fig.  5.33.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  Accord inq  to  Model  3  at  a  Distance  of 
1  km,  2  ..sec  Rise  Time,  and  25  ..sec  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fiq.  5.34.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  According  to  Model  3  at  a  Distance  of 
5  km,  2  usc-c  Rise  Time,  and  25  -sec  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fiq.  5.35.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  Accordina  to  Godel  3  at  a  Distance  of 
10  km,  2  ..sec  Rise  Time,  and  25  nsec  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fig.  5,36.  Computer  Solutions  to  Equation  (39), 
the  Electric  Field,  Accord im  to  Model  3  at  a  Distance  of 
100  km,  2  usee  Rise  Time,  and  2  5  ..sec  Fall  Time  for  the 
"'.round  and  at  Five  Altitudes. 
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Fiq.  5.37.  Computer  Solutions  to  Equation  (40), 
the  Maanetic  Field,  Accordinq  to  Model  1  at  a  Distance  of 
1  km,  500  sc c  Rise  Time,  and  10  ::sec  Fall  Time  tor  the 
Ground  and  at  Five  Altitudes. 


on  rt 


Ground 


8 , QOQ  feet 


Tig.  5.38.  Computer  Solutions  to  Equation  (40), 
he  Maonetic  Field,  /according  to  Model  1  at  a  Distance  of 
km,  500  .;sec  Rise  Time,  and  10  ,:sec  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fir:.  3.39.  Computer  Solutions  to  Equation  (40), 
ic  Field,  According  to  Model  1  at  a  Distance  o 
;.sc>c  Rise  Time,  and  10  ..sec  Fall  Time  for  the 
and  at  Five  Altitudes. 
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die:.  5.40.  Computer  Solutions 
he  Magnetic  Field,  Accordina  to  Model 
00  km,  500  ..sec  Rise  Time,  and  10  ..so 
round  and  at  Five  Altitudes. 

to  Eauation  (40) , 

1  at  a  Distance  o 
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r  u: .  5.42.  Computer  Solutions  to  Equation  (40), 
o  ianctic  Field,  Accord  inci  to  Model  1  at  a  Distance  o 
Km,  1  ..sec  Rise  Time,  and  25  ,;sec  Fall  Time  for  the 
ound  and  at  Five  Altitudes. 
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5.43.  Computer  Solutions  to  Equation  (40), 
Field,  Accordinq  to  Model  1  at  a  Distance  of 
Rise  Time,  and  25  ..sec  Fall  Time  for  the 
Five  Altitudes. 
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Fiq.  5.44.  Computer  Solutions  to  Equation  (40) 
netie  Field,  Accordinq  to  Model  1  at  a  Distance  ' 
1  ..sec  Rise  Time,  and  35  ..sec  Fall  Time  for  the 
and  at  Five  Altitudes. 
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Fig.  5.46.  Computer  Solutions  to  Equation  (40), 
;ho  Magnetic  Field,  According  to  Model  1  at  a  Distance  o 
i  Km,  2  ;,sec  Rise  Time,  and  50  usee  Fall  Time  for  the 
1  round  arid  at  Five  Altitudes. 
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Fic:.  5.47.  Computer  Solutions  to  Equation  (40), 
the  Magnetic  Field,  According  to  Model  1  at  a  Distance  o 
10  km,  2  ..sec  Rise  Time,  and  50  ;;sec  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fiq.  5.48.  Computer  Solutions  to  Equation  (40), 
the  Maq net ic  Field,  According  to  Model  1  at  a  Distance  of 
100  km,  2  usee  Rise  Time,  and  50  usee  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fig.  5.49.  Computer  Solutions  to  Equation  (40) 
the  Magnetic  Field,  Accordinq  to  Model  2  at  a  Distance  of 
1  km,  500  nsec  Rise  Time,  and  10  ,:sec  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fiq .  5.50.  Computer  Solutions  to  Equation  (40), 
he  Maanetic  Field,  Accordinq  to  Model  2  at  a  Distance  of 
km,  500  nsec  Rise  Time,  and  10  ..sec  Fall  Time  for  the 
'.round  and  at  Five  Altitudes. 
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Fiq.  5.51.  Computer  Solutions  to  Equation  (40), 
the  Maqnetic  Field,  Accordina  t1"1  Model  2  at  a  Distance  of 
10  km,  500  nsec  Rise  Time,  and  10  ..sec  Fall  Time  for  the 
Ground  and  at  Five  /Altitudes. 
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Fig.  5.52.  Computer  Solutions  to  Equation  (40), 
the  Magnetic  Field,  According  to  Model  2  at  a  Distance  of 
100  km,  500  nsec  Rise  Time,  and  10  usee  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fig.  5.53.  Computer  Solutions  to  Eauation  (40), 
the  Magnetic  Field,  According  to  Model  2  at  a  Distance  of 
1  km.  1  ..sec  Rise  Time,  and  25  ..sec  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fiq.  5.54.  Computer  Solutions  to  Equation  (40), 
the  Maonotic  Field,  Accordinq  to  Model  2  at  a  Distance  o 
5  km,  1  '..sec  Rise  Time,  and  25  ..sec  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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■  .‘lionet  ic  Field,  According  to  Model  2  at  a  Distance  o 
fm,  1  .  sec  Rise  Time,  and  25  ..sec  Fall  Time  for  the 
rand  and  at  Five  Altitudes. 
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Fici.  5.57.  Computer  Solutions  to  Equation  (40), 
the  Magnetic  Field,  According  to  Model  2  at  a  Distance  of 
1  km,  2  .isec  Rise  Time,  and  50  i:sec  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fiq.  3.58.  Computer  Solutions  to  Equation  (40), 
the  Magnetic  Field,  Accordina  to  Model  2  at  a  Distance  of 
5  km,  2  nsec  Rise  Time,  and  50  nsec  Fall  Time  for  the 
Ground  and  at  Five  /altitudes. 
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Fig.  5.59.  Computer  Solutions  to  Equation  (40) , 
the  Magnetic  Field,  According  to  Model  2  at  a  Distance  of 
10  km,  2  .  sec  Rise  Time,  and  50  ;.soc  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fia.  5.61.  Computer  Solutions  to  Equation  (40), 
ho  Mac  nr- tic  Field,  According  to  Model  3  at  a  Distance  o 
km,  500  nsec  Rise  Time,  and  10  -sec  Tall  Time  for  the 
round  and  at  Five  Altitudes. 
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Fig.  5.62.  Computer  Solutions  to  Equation  (40), 
the  Magnetic  Field,  According  to  Model  3  at  a  Distance  o 
5  km,  500  nsec  Rise  Time,  and  10  ,;sec  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fin.  5.63.  Computer  Solution?  to  Equation  (40), 
the  Maqnet ic  Field,  Accordina  to  Model  3  at  a  Distance  o 
10  kni,  500  nsec  Rise  Time-,  and  10  ..sec  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fiq.  5.64.  ConiDuter  Solutions  to  Eauation  (40), 
the  Magnetic  Field,  According  to  Model  3  at  a  Distance  of 
100  km,  500  nsec  Rise  Time,  and  10  a  sec  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fiq.  5.05.  Computer  Solutions  to  Equation  (40) 
the  Magnetic  Field,  Accordinq  to  Model  3  at  a  Distance  < 
1  km,  1  :.soc  Rise  Time,  and  25  usee  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fiq.  5.66.  Computer  Solutions  to  Equation  (40), 
the  Magnetic  Field,  According  to  Model  3  at  a  Distance  o 
5  km,  1  ..sec  Rise  Time,  and  25  ..see  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fin.  5. OR.  Computer  Solutions  to  Equation  (40), 
he  Marmot  ic  Field,  Accord  inc?  to  Model  3  at  a  Distance  o 
00  km,  1  ..sec  Rise  Time,  and  25  ..sec  Pal]  Time  for  the 
round  and  at  Five  Altitudes. 
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Piq.  5.09.  Computer  Solutions  to  Equation  (-4  0 
otic  Field,  Accord  inn  to  Model  3  at  a  Distance  of 
..sec  Rise  Time,  and  50  ,:sec  Tall  Time  for  the 
nu  at  Five  Altitudes. 
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Fir;.  5.70.  Computer  Solutions  to  Equation  (40) 
tno  Magnetic  field,  Accord inq  to  Model  3  at  a  Distance 
5  km,  2  .  sec  Rise  Time,  and  50  ..sec  Tall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fio.  5.71.  Computer  Solutions  to  E'.quation  (40), 
the  Magnetic  Field,  According  to  Model  3  at  a  Distance  o 
10  km,  2  usee  Rise  Time,  and  50  usee  Fall  Time  for  the 
Ground  and  at  Five  Altitudes. 
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Fig.  5.72.  Computer  Solutions  to  Equation  (40) , 
he  Magnetic  Field,  Accordinq  to  Model  3  at  a  Distance  of 
00  km,  2  usee  Rise  Time,  and  50  usee  Fall  Time  for  the 
round  and  at  Five  Altitudes. 


derivative  gains  to  zero  after  thirty  to  forty  microseconds 
therefore,  the  induction  term  tends  to  predominate  during 
the  impulse  period. 

Experimental 

In  comparing  the  theoretical  data  to  experimental 
data,  what  was  looked  for  was  similarity  in  the  general 
waveshape.  The  overall  magnitude  was  not  considered 
important  since  the  peak  current,  I  ,  is  a  multiplying  con¬ 
stant  that  can  be  pulled  out  of  the  integrals  of  Equations 
(39)  and  (40) .  This  means  that  if  a  model  predicts  a  field 
with  a  shape  closely  corresponding  to  experimental  data, 
but  whose  overall  magnitude  is  off,  the  predicted  field  can 
be  easily  amended  by  simply  increasing  or  decreasing  Iq 
accordingly.  The  value  of  IQ  needed  to  predict  both  the 
shape  and  magnitude  of  experimental  data  does  become  impor¬ 
tant  when  the  IQ  necessary  to  closely  approximate  data 
falls  outside  the  range  of  typical  values— ten  to  twenty 
kiloamps--as  is  the  case  in  the  second  model  (Lin  et  al., 
1980)  . 

Data  from  the  measurements  of  the  electric  and 
magnetic  fields  of  two  return  strokes  (not  of  the  same 
flash)  are  used  in  this  thesis.  The  data  for  the  first 
event  includes  simultaneous  electric  field  measurements  on 
the  ground  (Figure  5.73)  and  at  1500  feet  (Figure  5.74) 
at  a  distance  of  ten  kilometers  from  the  channel.  Data  for 
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the  second  event  consists  of  only  an  electric  field  measure¬ 
ment  at  8300  feet  (Figure  5.75) .  All  data  used  in  this 
thesis  was  measured  in  an  analog  FM  channel  with  500  kHz 
bandwidth  in  either  a  ground  station  or  an  Air  Force  WC-130 
aircraft  sponsored  by  the  Flight  Dynamics  Laboratory  at 
Wr ight-Patterson  AFB,  Ohio. 

Other  data  at  this  writing  is  still  being  pro¬ 
cessed.  Also,  other  data  was  processed  and  available  but 
did  not  closely  fit  the  points  at  which  the  fields  were 
theoretically  predicted.  Data  was  also  taken  with  a  VHF 
system;  this  data  was  preferred  over  the  FM  data  but  was 
unavailable  for  use  in  this  thesis. 

.  Looking  at  Figures  5.73,  5.74,  and  5.75  and  taking 
into  account  and  accepting  the  possibility  of  data  inver¬ 
sion,  the  modified  Lin  model  used  by  Master  et  al .  is 
clearly  the  best  fit.  However,  the  only  conclusion  that 
could  be  drawn  from  this  comparison  is  that  Lin's  modified 
model  as  presented  in  Master  et  al . ,  1981  best  fits  experi¬ 
mentally  measured  data  at  ten  kilometers  for  ground  level, 
rifteen  hundred  feet,  and  eight  thousand  feet;  the  question 
of  model  performance  for  other  distances  and  altitudes 
still  remains  unanswerable  at  this  time.  Even  at  that,  a 
conclusion  of  this  nature  is  not  well  supported  due  to  the 
limited  amount  of  experimental  data.  Clearly,  more  data 
for  comparison  purposes  is  needed  to  draw7  a  valid  conclusion 
of  the  performance  of  these  four  models. 


jioaei  1 


Model  2 


if 


Model  3  Master  et  al.,  1981 


Fig.  5.74.  Comparison  of  the  Electric  Field  Pre¬ 
dicted  by  the  Three  Models  and  Lin's  Modified  Model 
(Master  ct  al.f  1981)  with  Experimental  Data  Taken  at 
1500  Feet  at  a  Distance  of  10  kilometers.  Estimated  Rise 
Time  of  2  ;isec. 
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Fiq.  5.75.  Comparison  of  the  Electric  Field  Pre¬ 
dicted  by  the  Three  Models  and  Lin's  Modified  Model 
(Master  et  al.,  1981)  with  Experimental  Data  Taken  at 
8300  Feet  at  a  Distance  of  10  kilometers.  Estimated  Rise 
Time  of  2  :.;scc. 
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VI . 


Conclusions  and  Recommendations 


Concl us  1 ons 

The  only  conclusion  that  can  be  drawn  at  this  time 
is  that  the  modified  version  of  Lin's  model  (Master  et  al., 
1981)  best  predicts  the  measured  data  at  a  distance  of  ten 
kilometers  at  ground  level,  fifteen  hundred  feet,  and  eight 
thousand  feet.  For  other  distances  and  other  altitudes 
more  data  is  needed. 

Recommendations 

The  first  recommendation  is  to  obtain  more  of  the 
data  that  the  Flight  Dynamics  Laboratory  of  Wright- 
Patterson  AFB ,  Ohio  collected  in  Florida  in  August  1981. 

At  this  writing,  much  of  the  data  was  unavailable  for  use 
in  this  thesis  due  to  incomplete  processing,  not  fitting 
the  theoretically  predicted  data  points,  or  poor  resolu¬ 
tion.  Clearly,  more  data  is  needed  in  order  to  draw  a  good, 
valid  conclusion  as  to  which  model  performs  best  under 
various  conditions. 

Another  recommendation  is  to  analyze  the  three  cur¬ 
rents  of  Lin's  model  separately.  This  would  give  more  con¬ 
clusive  evidence  as  to  which  currents  predominate  at  vari¬ 
ous  times  in  the  calculations.  Along  those  same  lines,  it 
may  prove  fruitous  to  analyze  each  of  the  currents  separately 
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in  terms  of  the  electrostatic ,  induction,  and  radiation 
fields  they  produce. 

A  third  recommendation  is  to  theoretically  analyze 
these  models  at  distances  closer  than  one  kilometer  and 
at  distances  between  ten  and  one  hundred  kilometers.  This 
would  not  only  allow  the  use  of  more  of  the  experimental 
data  but  provide  more  insight  to  the  more  affecting  (to 
the  aircraft)  closer  lightning  flashes. 

Another  recommendation  is  to  do  some  analysis  of 
the  models  for  faster  rise  times.  Recent  data  has  shown 
that  some  return  strokes  may  have  rise  times  as  fast  as 
one  hundred  nanoseconds  (private  conversation  with  Dr. 

Pete  Rustan,  Fall  1981) . 

Final  recommendations  are  to  put  some  sort  of  decay 
or  attenuation  on  the  leader  current  in  Lin's  model  and  to 
model  the  channel  with  an  unmatched  load  at  the  top. 
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